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ABSTRACT 

We report the discovery of seven strongly lensed Lyman break galaxy (LBG) candidates at z ~ 7 
detected in Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) imaging of Abell 1703. 
The brightest candidate, called A1703-zDl, has an observed (lensed) magnitude of 24.0 AB (26<r) 
in the WFC3/IR F160W band, making it 0.2 magnitudes brighter than the zgso-dropout candidate 
recently reported behind the Bullet Cluster and 0.7 magnitudes brighter than the previously brightest 
known z ~ 7.6 galaxy, A1689-zDl. With a cluster magnification of ~ 9, this source has an intrinsic 
magnitude of Hiqq = 26 A AB, a strong zs5o — J\2b break of 1.7 magnitudes, and a photometric 
redshift of z ~ 6.7. Additionally, we find six other bright LBG candidates with Hiqq band magnitudes 
of 24.9 — 26.4, photometric redshifts z ~ 6.4 — 8.8, and magnifications \i ~ 3 — 40. Stellar population 
fits to the ACS, WFC3/IR, and Spitzer /IRAC data for A1703-zDl and A1703-zD4 yield stellar masses 
(0.7 — 3.0) x 10 9 M , stellar ages 5 — 180 Myr, and star- formation rates ~ 7.8 M yr _1 , and low 
reddening with Ay < 0.7. The source-plane reconstruction of the exceptionally bright candidate 
A1703-zDl exhibits an extended structure, spanning ~ 4 kpc in the z ~ 6.7 source plane, and shows 
three resolved star-forming knots of radius r ~ 0.4 kpc. 

Subject headings: galaxies: high-redshift — gravitational lensing: strong — galaxies: clusters: indi- 
vidual: Abell 1703 



1. INTRODUCTION 



The recently installed Wide Field Camera 3 (WFC3) 
aboard the Hubble Space Telescope has led to a signifi- 
cant increase i n the sample of z > 7 galaxy candidates in 
the past year (|Oesch et al.||2010b| [ Bouwen s et al.|20 10a; 
Bunker et al.|2010| IMcLure et al.|2UlU||Finkelsteiii et al. 



20101 |Bouwens et al.|2011b||Trenti et al.|2011||Yan et a 



2011). Already, more tha n 132 z ~ 7 — 8 Lyma n-break 



galaxy (LBG) candida tes (Bouwe ns et all2011b| see also 
Lorenzoni et al.||2010| [McLure et al. " 201l| ) have been 
found in ultradeep WEC3/1R observations of the Hub- 



ble Ultra-Deep Field (HUDF) and nearby fields, with 
even a few candidates at z ~ 8.5 and one at z ~ 10 
(Bouwens et al.|[2011a). These observations provide our 
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first glimpse of galaxies during the reionization epoch, 
showing a rapidly evolving galaxy luminosity function 
(LF) and a d eclining star- for mat ion r ate with increas- 



ing redshift (Bouwens et al 
of these z > 7 galaxies have 



2011b|a[ ). Recent studies 
also looked at their struc- 



ture and morphologies (Oesch et al. 2Q10a|, rest-frame 
i^-continu um slopes (Bouwens et al. 20105} |Finkelstein] 
| et al. | [20To| a nd star-formation rates and stellar masses 
( |Labbe et aL|2010b|a| |McLure et aT|2011 ). 

The ultradeep observations are complemented by shal- 
lower wide-field WFC3/IR surveys for z > 7 galaxies 
such as the Brig htest of Reionizing Galaxies (BoRG) 
flTrenti et aL]|2011| | and Hubble Infrared P ure Paral- 



lel Im aging Extragalactic Survey (HIPPIES) ( Yan et al. 
2011), which so far have uncovered 4 bright (25.5 — 26.7) 



z > 7.5 galaxies over 130 arcmin 2 . Over the next three 
years, the Cluste r Lensing And Super nova survey with 
Hubble (CLASH; postman et al.||2011[ ) and Cosmic As- 
sembly Near- infrared Deep Extragalactic L egacy Survey 
(CANDELS; IGrogin et al.l[20TT| |Koekemoer' etaTpQlT ) 
Multi-Cycle Treasury (MC I) programs will further aug- 
ment the sample of bright z > 7 galaxy candidates and 
our understanding of the high-redshift universe. These 
wide-field surveys are needed to characterize the bright 
end of the galaxy LF, where bright z > 7 galaxies are 
rare. Placing tighter constraints of the number density 
of bright sources also helps break the degeneracy between 
the characteristic luminosity L* and the faint-end slope 
a, a parameter whose value is crucial in determining the 
contribution of galaxies to the reionization of the uni- 
verse. 

The use of massive galaxy clusters as "cosmic" grav- 
itational telescopes has uncovered some t he brightest 
(< 26.5) z > 5 high-redshift galaxies to date (Franx et al. 
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1997; |Frye et al.|[2002] |Kneib et al.||2004; Egami et al. 



2005; Bradle y"etal. 1 120081 IZheng et al.| |2009; 
eraLlf2009[ |Hall et al.||2011|) and some of the 



Bouwens 
most dis- 



tant galaxies kno wn at the time of their discovery dFranx 



et al. 1997 iKneib et al.lpOl [Bradley et al. p008) 



Gravitational lensing by massive galaxy clusters can am- 
plify both the size and flux of background sources con- 
siderably. The increased spatial resolution allows high- 
redshift galaxies to be studied in unprecedented detail, 
provi ding clear views of their sizes and morphologies 
(e.g., |Franx et al.||1997l IKneib et al.||2004l [Bradley et al.| 
|2008 jzheng et al.|2009[]Swmbank et al.|2009p . This was 
clearly demonstrated with a very hi gh source-plane res- 
olution of 50 pc recently achieved by Zitrin et al. (2011) 
for the z = 4.92 galaxy behind MS 1358. Likewise, the in- 
creased brightness can place z > 7 galaxies within reach 
of ground-based spectroscopic follow-up observations. 

Luminous high-redshift galaxies are extremely valuable 
because their spectra can provide direct measurements 
of th e early star-form ation rate via Lya and Ha emis 



sion (lye et al. 2006) and the evolutio n of metallicity 
via metal e mission and absorption lines ( |Dow-Hygelund 
et al.|[2005 ). The spectra of z > 7 objects also pinpoint 



the epoch of the intergalactic medium (IGM) reionization 
through the effect of neutral h ydrogen in inhibiting the 
emission of Lya from galaxies ([Santos 2004; Malhotra 



Rhoads||2004| |Stark et al.||2010p . A truly neutraU^M 
will produce a damp ed Lya absorption profile (|Miralda- 



|Escude fc Rees|1998| ) that can be measured even at a low 
spectral resolution. 

Here we present the discovery of seven bright strongly 
lensed LBG candidates at z ~ 7 behind the massive 
galaxy cluster Abell 1703. The brightest candidate, 
A1703-zDl, is observed at 24.0 AB in the H 160 band, 
making it 0.2 magnitudes brighter than the zsso-dropout 
candidate recent ly reported behind the Bullet Cluster 
(Hall et al. 2011) and 0.7 magnitudes brighter than the 
previously brightest known z ~ 7.6 galaxy A1 689-zDl, 
found behin d the massive cluster Abell 1689 ( |Bradley| 
eTari|2QQ8| >. This paper is organized as follows! We 
present the observations and photometry in § 2 and 
dropout selection in § 3. In § 4 we discuss the source 
magnifications. We present the photometric redshifts in 
§ 5 and stellar population synthesis models in § 6. The 
results and the properties of the sources are discussed in 
§ 7. We summarize our results in in § 8. Throughout this 
work, we assume a cosmology with Vt rn = 0.3, = 0.7, 
and Hq = 70 km s _1 Mpc -1 . This provides an angular 
scale of 5.2 kpc (proper) arcsec -1 at z = 7.0. All m agni- 
tudes are expressed in the AB photometric system ( Oke 



1974). 



2. OBSERVATIONS AND PHOTOMETRY 

2.1. HST ACS and WFC3/IR Data 

We observed Abell 1703 [z = 0.284; |Allen eFl^[l992 ) 
with a single field of the Advanced Camera tor Surveys 
(ACS) Wide Field Camera (WFC) in 2004 November 
as part of an ACS GTO program to study five massive 
galaxy clusters (HST-GO10325). The observations cover 
a 3.4' x 3.4' field of view and consist of 20 orbits divided 
among six broadband filters: F435W (£^435; 7050 s), 
F475W (#475; 5564 s), F555W (V 555 ; 5564), F625W 
(r 62 5; 9834 s), F775W (i 775 ; 11128 s), and F850LP (z 850 ; 



17800 s). The ACS/WFC data were reduced wi th our 
ACS GTO APSIS pipeline flBlakeslee etlj][2003| >. The 
reductions reach 5cr limiting magnitudes (0.19" diameter 
aperture) of 28.5, 28.6, 28.2, 28.6, 28.4, and 28.0 in the 
£435, #475, ^555, ^625, ^775, and z 850 bands, respectively. 

We obtained WFC3/IR observations of Abell 1703 
in the F125W (J125) and F160W (H 160 ) bands, each 
with an exposure time of 2812 s, in 2010 April with 
the primary goal to search for z ~ 7 galaxies (HST- 
GO11802). The WFC3/IR observations cover the cen- 
tral 123" x 136" high-magnification region of the cluster 
(see Fig. [I). The depths of the WFC3/IR data reach 5<r 
limiting magnitudes of 27.3 and 26.9 in a 0.45" diameter 
aperture for the J125 and Hiqq bands, respectively. 

For the reduction of both the ACS and WFC3/IR data, 
we weight each individual exposure by its inverse vari- 
ance created from the sky background, modulated by the 
flatfield variations, along with the read noise and dark 
current. The drizzle combination procedure uses these 
inverse-variance images as weights and produces a final 
inverse- variance image for the combined, drizzled image 
in each filter. These inverse- variance weight i mages are 
used by SExtractor ( |Bertin fc Arnou t s 19 96) for both 
source detection and photometry (see §|2.3|). 



2.2. Spitzer/IRAC Data 

We utilized archival Spitzer/IRAC imaging of 
Abell 1703 (program 40311) obtained over three epochs 
between 2007 December and 2008 June. We used the 
Spitzer MOPEX calibration pipeline to combine the data 
in the 3.6 and 4.5 /xm bands over the three epochs. The 
total exposure times were 18.9 ks in the 3.6 and 4.5 /im 
bands, reaching 5<j limiting magnitudes of 24.7 and 24.1, 
respectively. 

2.3. Photometry 

We used SExtractor in dual-image mode for object 
detection and photometry. The detection image con- 
sisted of an inverse- variance weighted combination of the 
WFC3/IR J125 and H 160 images. We smoothed the ACS 
optical images to match the WFC3/IR images and mea- 
sured colo rs in small s calable Kron apertures (Kron fac- 
tor of 1.2; Kronp 980). We then correct the fluxes mea- 
sured in these smaller apertures to total magnitudes by 
using the flux measured in a larger Kron aperture (factor 
of 2.5) on the detection image. Likewise, we apply aper- 
ture corrections for light falling outside of the large Kron 
aperture using the tabulated encircled energies provided 
in the instrument handbooks. 

We were able to obtain Spitzer/IRAC fluxes for only 
two of the brighter and isolated sources. IRAC fluxes 
for these candida tes were obta i ned using the deblend- 
ing algorithm of [Labbe et al | ( |2006[ |2010b[ ). Briefly, 
this method involved using the higher resolution HST 
WFC3/IR images to create model IRAC images for the 
source and its nearby neighbors (assuming no differences 
between the structure or size of sources at 1.25 microns 
and IRAC wavelengths). We then vary the normaliza- 
tion of each model image, for both the source and its 
nearby neighbors, to fit the IRAC observations. Finally, 
we subtract the best-fit model profiles for the neighbors 
and perform photometry for the sources of interest in 
a 2.5"-diameter aperture. The IRAC errors include a 
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Figure 1. Color image (£85o-7i25^i6o) of the galaxy cluster Abell 1703 (z = 0.28). The locations of the high-redshift zsso-dropout 
candidate galaxies are marked by red circles and ellipses. The image field of view is 123 /x x 136 /x and is shown with P.A.= 152°. The white 
contours r epresent the critical curves at z ~ 7. The dashed cyan and green ellipses denote the regions where the strong lensing model of 
|Zitrin et al.| ( |2010| ) predicts counterimages for A1703-zD2 and A1703-zD5a/5b, respectively. 



component due to the modeling error. The modeling er- 
ror can be quite large for sources near the edges of the 
WFC3/IR image because there is no source template for 
IRAC sources found beyond the WFC3/IR field of view. 

3. SELECTION OF z ~ 7 z 850 -BAND DROPOUT 
CANDIDATES 

We search for z ~ 7 galaxies using a two-color zs5o~ 
dropout selection criteria, based on the magnitudes mea- 
sured in the small scalable apertures. We require candi- 
dates to have z^o — J125 > 0.7 and J125 — H160 < 0.5. 
In addition, they must be undetected (< 2a) in each op- 
tical ACS band, with not more than one band showing a 
> 1.5a detection. Further, candidates must be detected 



at > 5<j in the J125 band. In cases where an object is not 
detected in a particular band, we assign the object with 
the la detection limit to calculate object colors. 

Because our z 850 — J125 color criteria is slightly blue r 
than the ^850 — ^125 > 0.9 used by |Bouwens et al. (2011b), 
we effectively extend our redshilt selec tion "window to 
somewhat lower redshifts. Use of the IBouwens et al.l 



(2011b) color criteria, which was explicitly chosen to ex- 
clude source with redshifts z < 6.5, would eliminate only 
one of our candidates. 

Using this criteria, we identified seven zsso-dropout 
galaxy candidates with observed i?i60 band magnitudes 
between 24.0 — 26.4. The candidates are named in de- 
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A1703-ZD1 



A1703-ZD2 



A1703-ZD3 



A1703-ZD4 



A1703-zD5a 



A1703-zD5b 



A1703-zD6 



A1703-zD7 



Figure 2. Postage stamp cutout images of the high-redshift zsso-dropout candidate galaxies from the HST ACS and WFC3/IR data. 
The cutout images are 6 /x x 6", corresponding to 31.4 kpc on a a side at z = 7, and are shown with P.A.= 130°. As discussed in the text, 
A1703-zD5a/5b most likely represent two star-forming knots within a single source. 




creasing order of their brightness in the ifi6o band, with 
the exception of the close pair A1703-zD5a and A1703- 
zD5b. As discussed in the next section, A1703-zD5a/5b 
most likely represent two star-forming knots within a sin- 
gle source. We note that the A1703-zD5a component is 
detected in the V555 band at low significance (2.1cr), but 
not in the overlapping #475 and tq25 bands. Thus, we 
interpret the slight l/555-band detection as a likely statis- 
tical fluctuation. 

All of the candidates are clearly resolved with the ex- 
ceptions of our two faintest two candidates, A1703-zD6 
and A1703-zD7. A1703-zD6 is unresolved with SExtrac- 
tor stellarity parameter of 0.97. A1703-zD7 is somewhat 
extended, but only slightly resolved with a stellarity pa- 
rameter of 0.5. Because A1703-zD6 is unresolved, under 
normal circumstances it cannot be ruled out as a low- 
mass L,T dwarf star. However, this candidate has sub- 
sequently been confirmed to be at z = 7.045 with dee p 
Keck spectroscopic observations ( Schenker et al.pOll ). 

While the WFC3/IR data were taken alter the ACS 
optical data, we can rule we rule out supernovae as con- 
taminants to our sample because the sources are either 
resolved or in the case of the unresolved source A1703- 
zD6, spectroscopically confirmed at z — 7.045. 



The positions of these sources in the Abell 1703 data 
are shown in Figure [I] Their properties are listed in 
Table [I] and cutout stamps showing each of the sources 
are presented in Figure[2] The zsso - J125 and J125 — H\m 
colors of the candidates are illustrated in Figure [3] As 
seen in this figure, the zgso — J125 color of A1703-zD2 is 
exactly at our selection limit (0.7). While it could be in 
our z ~ 7 zgso-dropout sample as a result of photometric 
scatter, this candidate is completely undetected in the 
deep optical ACS data. Its somewhat bluer z^^o — J125 
color means it is at the lower-redshift edge of our ZS50- 
dropout selection window. This is completely consistent 
with its photometric redshift of z = 6.4 (see §|5|, making 
it our lowest redshift candidate. 

The best candidate, A1703-zDl, is an extremely bright 
zgso-dropout candidate, with a Hiqq magnitude of 24.0, 
that appears to be r esolved in three separate knots 
(see Fig. [2] and § 7.2). In Figure [4] we present a his- 
togram 01 both the observed and intrinsic (unlensed) 
Hiqo magnitudes compared with the 73 z ~ 7 candidates 
found in the HUDF09 and its two parallel field s and the 
WFC3/IR E arly Release Science observations ( Bo uwens] 
eTsIpoTTb ). 
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Table 1 

Observed Photometry of High-Redshift Candidates 



Candidate 



R.A. 



Dec 



^850 



Jl 



25 



GO 



3.6 |im 



4.5 /xm 



Zphot 



A1703-zDl 

A1703-zD2 

A1703-zD3 

A1703-zD4 

A1703-zD5a c ' d 

A1703-zD5b c 

A1703-zD6 e 

A1703-zD7 



13:14:59.4183 
13:15:06.5089 
13:14:58.3860 
13:15:07.1889 
13:15:07.7650 
13:15:07.7036 
13:15:01.0068 
13:15:01.2696 



51:50:00.843 
51:49:17.960 
51:49:57.740 
51:50:23.552 
51:49:09.333 
51:49:10.139 
51:50:04.353 
51:50:06.052 



25.8 ± 0.20 

25.6 ± 0.20 

26.8 ± 0.48 

> 28.0 

26.7 ± 0.34 
26.3 ± 0.29 

27.9 ± 0.53 

> 28.5 



24.1 ± 0.04 
24.9 ± 0.10 
25.5 ± 0.14 

25.5 ± 0.10 

25.6 ± 0.12 
25.3 ± 0.11 
25.8 ± 0.08 
26.8 ± 0.22 



24.0 ± 0.06 
24.9 ± 0.14 

25.1 ± 0.15 
25.4 ± 0.13 
25.7 ± 0.18 

25.3 ± 0.15 
25.9 ± 0.12 

26.4 ± 0.21 



23.9 ±0.1 24.7 ±0.4 



25.6 ± 0.5 24.7 ±0.5 



q + 0.9 

U -4.4 

o+43.6 

8 -10.4 

+ 1.3 

-2.5 

-1+0.2 

± -1.0 

n + 26.9 

U -15.0 

q+19.9 

^-8.0 

9 + 0.3 



6.7! 



7+0-2 
-0.1 

6 4+ - 1 

u, ^-0.3 



6.7 
8.4 
6.5 
6.5 
7.0 



+i.< 

0.2 
+ 0.9 
1.4 
+ 0.2 
0.3 
+ 0.2 
0.2 
+ 0.6 
0.2 
+ 0.2 
-1.7 



Note. — The sources without quoted IRAC magnitudes either suffer from significant confusion from neighboring sources or do not 
show an especially prominent (> 2a) detection. 

a The magnification errors represent the extreme values obtained from the minimum and maximum magnifications obtained within ±0.5" 
of each candidate and assuming Az ± 1.0 for the source redshifts. 

b Photometric redshifts determined from the BPZ code (Benitez 2000J. Because of the limited depth of the optical data, there is a small 
chance that some of the sources could be at low redshift (see §|o|r 

c As discussed in the text, A1703-zD5a/5b most likely represent two star-forming knots within a single source. The magnification of the 
combined A1703-zD5 source is 31.9l^;|. 

d This component of the zD5 candidate is detected in the V555 band at low significance (2.1cr), but not in the overlapping #475 and r625 

bands. Thus, the slight Vsss-band detection is likely a statistical fluctuation. 

e This candidate is spectroscopically confirmed to be at z = 7.045 Schen ker et al."| ( |2011| ). 




0.0 0.5 

(^125 — ^16())aB 

Figure 3. Z850 — J125 vs. J125 — -Hi 60 two-color diagram used 
to select the zsso-band dropout candidates. The error bars and 
lower limits are lcr. The gray region represents the Z850 — J125 and 
J125 — if 160 colors of the selection criteria. The blue, green, and 
red lines show the expected colors of star-forming galaxies with 
UV-continuum slopes of (3 = -1, -2, and -3 (f\ oc A^), respectively. 
The cyan, magenta, and yellow lines depict the expected colors 
of low-redshift interloper galaxies. The green region indic ates the 
colors of low-mass L,T dwarf stars (e.g., |Knapp et al.|2004| . While 
there is considerable overlap in the colors ol L,T stars, only one of 
the Z850-dropout candidates (A1703-zD6) is unresolved. However, 
this candida te has been s ubsequently spectroscopically confirmed 
at z = 7.045 |Schenker et al.| ( |2011| >. 

4. SOURCE MAGNIFICATIONS AND COUNTERIMAGES 

Several detailed studies to model the lensing of 
Abell 1703 have been performed in recent years 



Lim ousin et al.||2008| |Richard et al. 2009 
|2U!Up . We adopt the |Zitrin et aTf 



Zitrin et al. 
Abell 1703 



strong lensing model to estimate the magnifications of 
the seven z ~ 7 source s and to identify possible coun- 
terimages. Zitrin et al. (2010) used 16 multiply-imaged 
systems behind Abell 1703 and applied two independent 
strong lensing techniques to the high-quality, multiband 
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Figure 4. Histogram of the observed and intrinsic (unlensed) 
if 160 magnitudes for the seven zg 50 -dropout [z ~ 7) candidate 
galaxies identified behind Abell 1703 (green and red) compared 
with the 73 z ~ 7 candidates found in the HUDF09 and its two 
parallel field s and the WFC3/IR E arly Release Science observa- 
tions (blue) jBouwens et al.|2011b| >. 

ACS data, yielding similar results. Their strong lensing 
model places tight constraints on the inner mass pro- 
file, and thus provides reliable magnification estimates 
for background sources. The magnifications of the high- 
redshift candidates range from \i ~ 3 to large magnifica- 
tions of ~ 25 — 40, found for three of our sources that are 
located near the critical curve, where the magnification 
formally diverges. The magnification of each candidate 
is presented in Table [l] 

We estimated the magnification uncertainties by tak- 
ing models extracted from the la confidence level, as de- 
termined by the x 2 minimization of model, and marginal- 
izing over the true la errors. To make the error estimates 
more conservative, we also incorporate the range of mag- 
nifications obtained within ±0.5" of each candidate and 
apply a Az ±1.0 to the redshift of each source. Thus, 
the ±0.5" shift is a measure of the magnification variance 
around the location of the source and the application of 
Az ±1.0 accounts for the possible local uncertainty in 
the location of the critical curves. For objects that are 
close to the critical lines, the magnification errors are di- 
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verging due to their proximity to the critical curve, while 
objects far away will have a well-determined magnifica- 
tion as the latter slowly varies in regions away from the 
critical curve. 

The brightest candidate, A1703-zDl, has a magnifica- 
tion of ~ 9, giving it an intrinsic magnitude of ~ 26.4 in 
the Hiqo band. The Abell 1703 strong lensing model pre- 
dicts counterimages for the three high-magnification can- 
didates, A1703-zD2 and the A1703-zD5a/5b pair, which 
are located nearby or on the high-redshift critical curve 
(see Fig. [T]). The lensing model predicts three counter- 
images for A1703-zD2 (/i = 24.8; see Fig.[l]). Taking into 
account the much smaller magnifications (/i = 5.5 — 9.0) 
of the counterimages, they are predicted to have Hiqo 
magnitudes between 26.0 — 26.5. This is sufficiently 
bright that there was some possibility that we might lo- 
cate them, but also a good chance we might not because 
they could easily be lost in the wings of a foreground 
galaxy. Despite an extensive search, we did not find any 
viable z ~ 7 candidates near the predicted positions of 
the counterimages. 

The close pair of Al703-zD5a and Al703-zD5b are also 
located in very close proximity of the critical curve and 
as such have high magnifications of \i ~ 27 — 40. These 
candidates are also predicted to have counterimages on 
the other side of the brightest cluster galaxy (BCG; see 
Fig. [I]) with magnifications of ji = 5.5, about 5 times 
less than A1703-zD5b. The predicted counterimages are 
expected to have an H\qo magnitude of > 27.0, which is 
fainter than our 5<r limiting magnitude of 26.9. Hence, 
it is not surprising that no ^850-dropout candidates are 
found in the predicted region of the counterimages. 

The critical curve lies only 2.6" from passing between 
the A1703-zD5a/b sources, which is within the model 
(and redshift) uncertainty. To test the hypothesis that 
these two candidates are multiple images of the same 
source, we constructed a new model assuming that zD5a 
and zD5b are the same object. The resulting model is 
physically plausible and the predicted counterimage of 
this system, located in the same region marked in Fig. [TJ 
is again fainter than the 5cr limiting magnitude of 26.9. 
Because the overall reproduction of all other systems re- 
mains the same, we cannot exclude this option based 
solely on the mass model. This possiblity could also be 
supported by their similar colors, morphologies, and pho- 
tometric redshifts. 

However, because surface brightness is conserved in 
lensing, zD5a and zD5b should have the same surface 
brightness if they are indeed the same source. For zD5a 
and zD5b, we find surface brightnesses of 25.9 and 25.8 
mag arcsec -2 , respectively, in the J125 band and 26.4 
and 25.9 mag arcsec -2 , respectively in the Hiqq band. 
We note tha t these results are consistent with Oesch 
al. (2010a) who found a mean Ji25-band observed 
surface brightness of ~ 26 mag arcsec -2 for a sample 
of z ~ 7 zgso-band dropouts candidates spanning 26 to 
29 mag in the J125 band. While their surface bright- 
ness agrees in the J125 band, zD5b has a brighter sur- 
face brightness in the Hiqo band. On these grounds we 
conclude that zD5a and zD5b are two unique sources. 
Further, their very close proximity of 0.76" in the image 
plane, with a magnification of 8.6 along the line sepa- 
rating them, translates to only ~ 480 pc in the source 



plane. Thus, even though we do not observe a diffuse 
component, which could have very low surface brightness 
below our detection limits, between them, we conclude 
that zD5a and zD5b are most likely two star-forming 
knots within a single source. The magnification of the 
combined A1703-zD5 source is 31.9t 2 o j. We refer to the 
these sources separately because they appear as distinct 
sources in our catalog. The small elliptical Kron aper- 
tures used to measure their colors are well separated with 
sizes of 0.34" x 0.26" and 0.28" x 0.22", respectively. 

We also considered the possibilty that zD5a/b and zD2 
were all counterimages of the same source. With zD5a/b 
representing two bright star-forming knots in a single 
candidate galaxy and not counterimages of the same 
source, we conclude that zD2 cannot be a counterim- 
age of the zD5a/b source based simply on morphology. 
This was verified by constructing an additional model 
considering the center of zD2 and the center of zD5a/b 
as counter positions of the same source. 
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Figure 5. Probability distributions of the photometric redshifts 
for each of the candidates. The vertical yellow line represents the 
redshift of the Abell 1703 cluster (z = 0.28). As discussed in the 
text, A1703-zD5a/5b most likely represent two star-forming knots 
within a single source. The A1703-zD5a component shows the 
highest probability of being at low redshift due to its modest 2.1a 
detection in the F555W band, which we attribute to a statistical 
fluctuation as described in the text. 



5. PHOTOMETRIC REDSHIFTS 

To estimate the redshifts of the candidates, we used 
the Bayesian photometric redshift (BPZ) code ([Bemtez 



2QQQ| |Bemtez et al.|[2004| |Coe et al.|[2QQ6). Briefly, the 



photometric redshifts are based on a \ l fitting procedure 
to template spectra. Because the shape of the redshift 
distribution is not well calibrated at z ~ 7, we assumed 
a flat prior for all redshifts. The photometric redshifts of 
the zgso-dropout candidates are presented in Table [l] and 
their posterior P(z) probability distributions are shown 
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in Figure |5j We find redshifts in the range of z p hot — 
6.4 — 8.8, with a median redshift of 6.7. 

The posterior probability distributions show that there 
is a small chance that A1703-zD3, A1703-zD5a, and 
A1703-zD7 could be at low redshift. Because of the 
modest 2.1a detection in the F555W band, the A1703- 
zD5a knot shows the highest probability of being at low 
redshift. While we cannot completely exclude the low- 
redshift solutions for these candidates, objects in this 
magnitude range (24.0—26.4) with colors similar to LBGs 
would be rare and are more likely to be at high redshift. 
The BPZ results indicate that the exceptionally bright 
candidate A1703-zDl has a narrow probability distribu- 
tion at z = 6.7 and has the highest probability of being 
at high redshift, not showing any evidence to suggest a 
low-redshift solution. 

6. STELLAR POPULATION MODELS 

We performed fits to the multiband HST and Spitzer 
photometry of A17Q3-zDl and A17Q3-zD4 us i ng the stel- 
lar populatio n mode l s of Bruzual fc Chariot (2003). We 
adopted a |Salpeter| ( 1955p initial mass function (IMF) 
with mass cutoffs of 0.1 and 100 M and models with 
subsolar (Z = 0.2Z© = 0.004) metallicities. The effect 
of dust reddening i s included in the models using the 



Calzetti et ah] |2000 ) obscuration law. We use thelMadaul 
( |l995| procedure to correct the models for Lyman-series 
line- blanketing and photoelectric absorption. The stellar 
population models are constrained such that the stellar 
age must be less than the age of the universe at the fit 
redshift (e.g., 0.75 Gyr at z = 7.0). We consider two 
star- formation histories (SFH): simple (single-burst) stel- 
lar population (SSP) models and constant star-formation 
rate (CSFR) models. 

The best-fit stellar population models for A1703-zDl 
and A1703-zD4, the two sources for which we were able 
to obtain IRAC photometry, are shown in Figures [6] and 
fn respectively and the parameters are given in Table § 
For these sources we find reasonably good model fits to 
the observed broadband photometry. For A1703-zDl, we 
note that the SED models are unable to fit the low flux 
in the IRAC 4.5/im band, resulting in a somewhat higher 
1.3-2.4 

We find intrinsic (unlensed) stellar masses for both 
candidates in the range from (0.7 — 3.0) x 10 9 M with 
star- formation rates of 7.3 ± 0.3 yr _1 and 8.2 ± 
1.2 M c — - 1 
z ~ 6 - 



yr % broadly consist ent with those found for 
alaxy candidates (Labbe et al. 2010b|a[ |Gon-| 



zalez et al. 2010; McLure et al.|201l[ ). Because these two 
candidates are located near the edge of the WFC3/IR 
field, there are no source templates for IRAC sources 
found beyond the WFC3/IR field of view. This results in 
a large modeling error for the neighboring sources, which 
we include in the total IRAC error for these candidates. 
While our modeling procedure weights each photometric 
data point by its inverse variance, because of the rela- 
tively large IRAC errors for these sources it should be 
noted that there is a larger uncertainty in the resulting 
stellar masses and ages derived from the stellar popula- 
tion models. 

W hile we assumed a Salpeter IMF, fitting models us- 
ing a |Chabrier (2003) IMF would yield lower masses and 
star- format ion rates by a factor ~ 1.5. For both candi- 
dates, we note the the CSFR models provide consistently 



older SFR-weighted mean stellar ages than that for the 
SSP models. We find CSFR models with weighted ages 
of 100 — 180 Myr, values again similar to that r ecently 
reported for a sample of z ~ 7 — 8 candidates ( Labbe 
etaLl|2010b| |McLure et aT^|20TT| ). 
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Figure 6. Best-fit stellar population models to the observed 
multiban d pho tometry of A1703-zDl. The models assume a 
|Salpeter| fl955| ) IMF with subsolar (Z = 0.2Z©) metallicity. The 
2(7 upper limits are shown for the ACS optical non-detections, but 
the models were fit to the measured fluxes and errors. The stellar 
masses are intrinsic values, corrected for the cluster magnification. 
The larger xi values (1.3 — 2.4) for these models results from the 
poor fit to the IRAC 4.5/mi band data. 
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Figure 7. Best-fit stellar population models to the observed 
multiband pho tometry of A1703-zD4. The models assume a 
|Salpeter| |l955| > IMF with subsolar (Z = 0.2Z©) metallicity. The 
2cr upper limits are shown for the ACS optical non-detections, but 
the models were fit to the measured fluxes and errors. The stellar 
masses are intrinsic values, corrected for the cluster magnification. 



7. DISCUSSION 



7.1. 



Al703-zDl Brightness 

Some of the brightest (< 26.5) z > 5 high-redshift 
1-19971 |Frye et al. 120021 [Kneib 



galaxies to date flFranx et ah] 



[etaT. 2004; E gami et al.|2UUk||Bradley et al.i2lJUgf [Z heng 
|et al.||2009| |Bouwens et al ||2009l |Hall et al.||2011| ) have 
been identified in searches behind strong lensmg clusters. 
Several of the more notable examples inclu de the bright 



^775-d ropout galaxy, A1703-iDl, found by Zh eng et al. 
|2009[ ) behind Abell 1703. This candidate has a N1C- 
MOS/NIC3 #i6o band magnitude of 23.9 and is lensed 
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Table 2 

Best-fit Stellar Population Model Results 



Candidate 


SFH a 


Zphot 


Mass c 
(1O 9 M ) 


Age w d 
(Myr) 


SFR 
(M© yr- 1 ) 


A v 


xl 


A1703-zDl 


SSP 
CSFR 


6.7 ±0.1 
6.86 ±0.1 


0.7 ±0.1 
1.5 ±0.05 


5.0 
101.3 


7.3 ±0.3 


0.7 ±0.1 
0.00 ±0.00 


1.3 
2.4 


A1703-zD4 


SSP 
CSFR 


8.4 ±0.3 
8.4 ±0.3 


2.0 ± 1.4 
3.0 ±0.4 


32.0 
180.1 


8.2 ±1.2 


09 +0 - 29 
0.0 ±0.0 


0.8 
0.8 



Note. — Models assume a |Salpeter| ( |1955| ) IMF with mass cutoffs of 0.1 and 100 M© 
and subsolar [Z — 0.2 Z@ = 0.UU4) metallicities. 

a Star-formation history: simple (single-burst) stellar population (SSP) or constant star- 
formation rate (CSFR) models. 

b Photometric redshifts derived from the stellar population fitting. These redshifts are 
consistent with those presented in Table [I] 

c Best-fit stellar mass, corrected by the magnification at the fitted Zph t redshift. 

d SFR- weighted mean stellar age (cf. |Forster Schreiber et al.||2004| . For a CSFR model, 

the weighted age is simply half of the time elapsed since the start of star formation. 



by a factor \i ~ 3.1. From the SED fitting of this source, 
we found a photometric redshift of z — 5.95 ± 0.15, 
which is consistent with t he Keck spectroscopi c redshift 
of z = 5.827 measured by [ Richard et al.| ([20 09 ) . 



At somewhat higher redshift, Kneib et al.|(|2004 ) iden- 
tified an exceptionally bright candidate (MC3 Hiqq= 
24.1) behind the galaxy cluster Abell 2218. This triply- 
imaged candidate, A2218-iDl, lies near the high-redshift 
critical curve and has a large magnification of /i ~ 25. 
While this very bright candidate has a similar i^i6o-band 
magnitude as A1703-zDl and has been suggested to be at 
z ~ 7, it has a zsso ~ Jno color of just ~0.4 magnitudes, 
which is more consistent with a redshift of z = 6.3 ±0.1. 
It therefore almost certainly has a lower redshift than 
the present candidate Al703-zDl and A1689-zDl, which 
is the brightest z ~ 7.6 candidate known. A1689-zDl 



was discovered by Bradley et al. ( |2008| ) behind the mas- 
sive galaxy cluster Abell 1689. This source is magnified 
by a factor of \i ~ 9.3 and has a NIC3 Hiqq magni- 
tude of 24.7, which is 0.7 magnitudes fainter than the 
exceptionally bright candidate z ~ 7 candidate A1703- 
zDl presented here. Most recently, |Hall et aL (2011) 
reported the discovery of a very bright Zgso-dropout can- 
didate behind the Bullet Cluster. The lensing model of 
the Bullet Cluster suggests that this candidate is doubly- 
imaged with observed i^i6o-band magnitudes of 24.2 and 
25.0 and magnifications of 8.4 and 12, respectively. How- 
ever, because this candidate is detected in the 2775 band, 
there is a f air chance this o bject could be a low-redshift 
interloper ( |Hall et al.|[2011[ ). 



7.2. AllOS-zDl Source-Plane Reconstruction and 
Morphology 

The large magnification of A1703-zDl allows us an op- 
portunity to examine the morphology of this exception- 
ally bright z ~ 6.7 galaxy candidate at very high spatial 
resolution. With a magnification of \i ~ 9, the strong 
lensing effect provides an increased spatial resolution by 
about a factor ~ 3 compared to an unlensed source. This 
permits us to resolve spatial structures that would oth- 
erwise be unobs ervable in high-reds hift z ~ 7 galaxies. 

1703 cluster 
the 



120101) Abell 
^03-zDl 



m 



source 



We used 
lensing model to reconstruct 
plane at z ~ 6.7. The deprojected image of A1703-zDl in 
the WFC3/IR J125 band is shown in Figure [8] The linear 




Figure 8. Source-plane reconstruction of A1703-zDl in the 
WFC3/IR J125 band. The source has an extended morphology that 
spans ~ 0.74" (~ 4 kpc) in the source plane at z = 6.7. A1703- 
zDl is comprised of three separate resolved knots, each with radius 
r ~ 0.08" (~ 0.4 kpc). Of course the physical size and structure 
we infer for this candidate is somewhat dependent on the details 
of the gravitational lensing model. The white ellipse denotes the 
WFC3/IR J125 PSF in the source plane (0.16" FWHM in the im- 
age plane). 

magnification along the shear direction is 4. 88^0' 14 an d 
1.84 ± 0.01 in the direction perpendicular to the shear 
direction. This candidate appears to be comprised of 
three resolved star-forming knots, each with a radius r ~ 
0.08" (0.4 kpc) in the source plane. Altogether, A1703- 
zDl has an extended linear morphology that spans ~ 
0.74" (~ 4 kpc) in the source plane at z = 6.7. Of course 
the physical size and structure we infer for this candidate 
is somewhat dependent on the details of the gravitational 
lensing model. 

There are now several examples of z > 5 lensed galaxy 
candidates that have morphologies consisting of star- 
forming knots ([Franx et aL 1997 Kneib et al. 2004 
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extended and show significant substructure. In partic- 
ular, both A1703-iDl and A1689-zDl show a pair of 
resolved star- forming knots. Additionally the pair of 
z ~ 6.5 dropout candidates, CL0024-iDl and CL0024- 
zDl, each seem to consist of two components. With a 
separation of only 2 kpc in the s ource plane and n early 
identical redshifts and properties (Zheng et al. 2009), it is 
possible that CL0024-iDl and CL0024-zDl are spatially 
associa ted or merging galaxies. Recently, |Oesch et al. 



d2010a[ ) even found extended features with resolved dou- 
ble cores in two unlensed z ~ 7 galaxies identified in 
the WFC3/IR HUDF. The apparent frequency of high- 
redshift galaxies showing substructure and multiple com- 
ponents is perhaps not surprising, but it provides strong 
evidence that both clumpy star formation and merging 
are important aspects of galaxy buildup at these very 
early epochs in the universe. The existence of substan- 
tial substructure is also expected based on studies of 



low-redshift Lyman-break analog galaxies (|Overzier et al. 
2008). 



7.3. Number Counts of z ~ 7 Candidates 

The discovery of seven z ~ 7 zgso-dropout candidates 
in a single cluster field is quite remarkable, especially 
since four of our candidates (zDl, zD3, zD6, and zD7) 
lie in a small region near the edge of the WFC3/IR im- 
age. While this may be surprising and call into question 
the reliability of these candidates, our confidence is bol- 
stered by the spectro scopic confirmation at z = 7.045 
flSchenker et al.|2011| of A1703-zD6, the second-faintest 
candidate. We consider the extremely bright (i^i60 = 
24.0) candidate A1703-zDl, with a strong z 850 — J125 
break of 1.7, to be a rather robust high-redshift candi- 
date, but we acknowledge that A1703-zD7, the faintest 
and nearly unresolved candidate, is certainly less secure 
than the other candidates. 

While the cluster magnification increases the effective 
depth of the observations, the source-plane area that is 
surveyed at high redshift decreases inversely proportional 
to the magnification factor. The Abell 1703 WFC3/IR 
image field of view is 4.6 arcmin 2 , but we estimate that 
we are effectively surveying only ~ 0.9 arcmin 2 in the 
z ~ 7 source plane. Thus, we derive a simple estimate 
of the number density of z ~ 7 sources in Abell 1703 
of 7.8 arcmin -2 . Typically, blank field surveys such as 



the HUDF09 and its two parallel fields (Bouwens et al. 
2011a| |Oesch et al.]2010b| ) find - 3.5 - 4.3 z~l sources 
per arcmin 2 . While our number density is larger than 
that found in typical blank fields, if we allow for the pos- 
sibility that zD7 is a low-redshift interloper, we are left 
with six sources or 6.7 arcmin -2 . However, we note that 
cosmic variance in the number counts of high-redshift 
source s is significant in these r elatively small- area HST 
fields (Trenti & Stiavelli |2008 ), especially the cluster 
fields with their significantly reduced area in the high- 
redshift source plane. 

We can place the observations of Abell 1703 in context 
with those obtained for a small, but growing, sample of 
strong lensing clusters with high-quality optical and NIR 
multiband data. While Hall et al. ( 2011[ ) has reported the 
discovery of ten 50 -drop outs behind the Bullet Cluster 
(Hall et al. 2011), most cluster lensing fields have pro- 
duced at most 1 — 2 z > 7 candidates. This includes 



NICMOS and WFC3/IR imaging of Abell 1689 and the 
recent 16-band imaging of Abell 383 and MACS1149 ob- 
tained by the CLASH MCT program. The apparent 
large variations in the number of z > 7 candidates dis- 
covered behind lensing clusters suggests, not surprisingly, 
that z > 7 galaxies may be highly clustered. One may 
therefore need to survey a large number of clusters to 
overcome the substantial large-scale structure effects. 

8. SUMMARY 

We report the discovery of a very bright, highly mag- 
nified LBG candidate (A1703-zDl) at z ~ 6.7 behind 
the massive galaxy cluster Abell 1703. A1703-zDl is 
0.2 magnitudes brighter than the recently disco vered 



£85o-dropo ut candidate behind the Bullet Cluster (Hall 
et al.|20lT| ) and 0.7 magnitudes brighter than the current 
brightest km 



lown z 



7.6 galaxy A1689-zDl, identified 
behind the massive galaxy cluster Abell 1689 ( |Bradley| 



et al 
least 



2008). We find a strong zgso — ^125 break ol at 
7 mag and best-fit photometric redshift of z = 6.7. 



Using the |Zitrin et al.| (|2010j) cluster lensing model, we 
estimate a magnitication oi /i = 9 at z ~ 6.7 at the posi- 
tion of A1703-zDl. The candidate is extended, spanning 
~ 4 kpc in the reconstructed source plane, and is resolved 
into three resolved star- forming knots. The source plane 
deprojection shows that the star formation is occurring 
in compact knots of size ~ 0.4 kpc. 

Additionally, we find six other bright z ~ 7 zsso- 
dropout galaxy candidates behind Abell 1703. One 
of these candidates, A1703-zD6, has been subsequently 
confirmed with Keck spectroscopy to be at z = 7.045 
Schenker et al.| ( |2011[ ) . The candidates are observed with 
Hiqq band magnitudes of 24.9 — 26 .4, with a wide range of 
magnifications from 3 — 40. Their photometric redshifts 
are found to be in the the range of z p hot = 6.4—8.8, with a 
median redshift of 6.7. Using stellar population models 
to fit the rest-frame UV and optical fluxes for A1703- 
zDl and A1703-zD4, we derive best-fit values for stellar 
masses (0.7 - 3.0) x 10 9 M , stellar ages 5 - 180 Myr, 
and star-formation rates ~ 7.8 M yr -1 . 

A1703-zDl, with a photometric redshift of z ~ 6.7 
is the brightest observed z ~ 7 galaxy candidate found 
to date. We are planning to observe A1703-zDl with 
near-IR spectroscopy to confirm its redshift and study 
its spectrum. Bright high-redshift galaxies such as these 
are valuable targets for ground-based spectroscopy and 
are pathfinding objects for future facilities such as the 
James Webb Space Telescope. 
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